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Introduction 


The following pages describe the results of work performed under NASA 
Consortium NCC2-5025, including the development of a program for structural 
analysis and weight estimation for an oblique all-wing aircraft. 

Any aircraft preliminary design study requires a structural model of the proposed 
configuration. The model must be capable of estimating the structural weight of a 
given configuration, and of predicting the deflections which will result from 
foreseen flight and ground loads. The present work develops such a model for the 
proposed Oblique All Wing airplane. The model is based on preliminary structural 
work done by Jack Williams and Peter Rudolph at Mdng, and is encoded in a 
FORTRAN program. As a stand-alone application, the program can calculate the 
weight CG location, and several types of structural deflections; used in conjuncbon 
with an aerodynamics model, the program can be used for mission analysis or sizing 
studies. 



Program Structure 

The program consists of a group of analysis and drawing subroutines which 
communicate with each other and with a database through the interface package, GENIE 1 . 
The database contains all physical parameters necessary to describe the airplane: element 
dimensions, material parameters, component locations, weights of fixed equipment, etc. 
Figure 1 shows the structure of the interaction between the program's subroutines and the 
database. 

Analysis routines 

• Plate BeamProp 1 calculates structural properties of plate-like and beam-like elements 

used to build up the wing section structure. 

• CompSec computes, for a given wing section, such properties as structural weight per 

unit span, section CG location (fore and aft), and wing stiffness (El). 

• Buildwing assembles section properties from CompSec into a single beam model for 

the whole airplane, adds in the user-defined fixed equipment weight distribution, 
and computes the total airplane weight, weight distribution, CG location, and 
bending stiffness. 

• WingSupport finds the forces to which the wing will be subjected, both in flight and 

when resting on its landing gear. 

• InternPress describes the cabin pressure, for use in calculation of the skin 

pressurization pillowing. 

• Deflection takes the results of all the above routines, and calculates the wing bending 

moments and deflections across the span. 

Drawing Routines 

• Planform draws the planform of the wing described in the database, including engine 

and fin locations; lines on the drawing marie locations where the wing section is 
explicitly defined in the database . 


’ Available from Desktop Aeronautics, PO Box 9937, Stanford, CA 94305. 

-Note that throughout the rest of this text, words in italic type face denote program subroutines or names in 
the database. 
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• InspectSection is an interactive routine which allows die user to graphically examine 

the parameters of the wing section he has defined in the database. 

• Sending produces a front view of the airplane, tncluding the bending deflections 

calculated by the analysis routines. 

Wing Description 

" K><l ' led “ 3 Slen<fcr b “ n ’ of sk ’ wl ) r ^8 cross section. The 
base defines NS, mom spanwtse stations (SuSuuionf), indexed from 1 to NSuuions 

from the left to nght across the span), where the moments and deflection are to be 
calculated. Structural properties are calculated at several (fewer than NSuuions) locations 
Where the wtng section is explicitly defined (see Section Description below), and linearly 
interpolated to the NSuuions stations. All wing properties are assumed to vary linearly 
between these stations, but may jump discontinuously across each station, requiring that 
two values be defined for each property at each station, a T and a V value (see figure 2) 
The following properties are calculated or interpolated at each station: 


StalElneniaO, StarElmeria() 
StaJLinDensO, StarUnDensO 
StalCGxQ, StarCGxf) 

StalCentroidxO, StarCentroidxO 
StalChondQ, StarChordO 
Srallexf), Starlexf) 

StahovercO, StartovercQ 
StalflapcQ, StarftapcO 
StalmcPillCoeffO, StarmcPillCoeffO 
StalmcBuckCoeffO, StarmcBuckCoeffO 


Wing stiffness distribution 
Weight distribution 
Chordwise local CG location 
Chordwise section centroid location 
Local wing chord 
Local leading edge position 
Local thickness to chord ratio 
Local flap chord 
Pillowing coefficient 
Buckling coefficient 


Section Description 

.. . f ch is dcfin « l ^ applying 48 parameters: 14 geometric parameters 

The u <SCn * ** airf ° l1 ' a ” d 34 ma *™ 1 Parame,CTS for *ctron components. 

The 14 geometric parameters define a rough outline of tire airfoil by designating tire 

chordwise location, thickness, and clramber (in unto of local chord, SecCuZ a. four 
stations (see Figure 3): 


SecfixovercO, Secfitovercf), SecfeCamber() 
SecCrestxovercO, SecCresttovercf), SecCrestCamberf) 
SecrsxovercO, Secrstovercf), SecrsCamberQ 
Secchxovercf), Secchtoverc (), SecchCamberQ 


Forward spar 
Crest 
Rear spar 
Control hinge 
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Any Wing Property 
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View: 



5S2? 2 the T ng Properties are explicitly defined or calculated at a number of 
require that a teftand a right^M^ * ff0 * )ert * es 
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The section’s location fore-and-aft is defined by its leading edge position in inches 
relative to the aircraft datum. 


Hex) 


The 34 material parameters describe two types of section elements (see Figure 4)- 

plate-like members (e.g. skins, spars, ribs), and beam-like members (e.g. spar-sldn 

spltces, nb-chords). Each plate-Uke element is defined by its outer thickness, skin 
thickness, and core type: 


l ec leShnOT(), SecleSkinSTO, SecleSkinCTf) 
^t^c S J^^(K SeancSkinST ^ SecmcSkinCTO 
\f l S J^? T J } ^l^ eStinST() ’ Se *eSkinCT() 
S c e< t9l a Se & ST 0, SecfiCTO 
SecfiOTQ, SecfsSTf), SecftCT() 

S c eC ?i R L^()’J ec ™ RibST 0. SecmcRibCTO 
SecteRibOTO, SeaeRibSTQ, SeaeRibCTQ 


Leading edge skins 
Midchord skins 
Trailing edge skins 
Forward spar material 
Rear spar material 
Midchord ribs 
Trailing edge ribs 


Tire outer thickness is defined as the total plate thickness, including all core material and 
Ming sheets. The skin thickness is taken to mean the thickness of the high density lacing 
sheet material attached to the outside of the core material. The plate core is either of macro- 
core or a titanium core material. 

Each beam-like element is described by a beam depth and an area ftaction: 


SecfiCapDO, SecfiCapAFQ 
SecrsCapDO, SecrsCapAFQ 
SecleRibTubeDO, SecleRibTUbeAFO 
SecleRibChordDO, SecleRibChordAFO 
SecmcRibChordAFf) ' 

SecteRibChordAFO 


forward spar- skin splices 
rear spar-skin splices 
leading edge rib tubes 
lead edge rib chords 
midchord rib-skin splices 
trailing edge rib chord 


Ahbeam elements are assumed to be made of resin transfer molding (RTM). The bending 

inertia of die wing section depends strongly on the amount of material in these beam-like 
elements, but only weakly on the beam's own bending inertia. The beam's enact geometry 
is therefore unimporian. for me present analysis, and the beam cross section can simply be 
characterized by the beam size and material ftaction. 

Three parameters describe the rib spacing in inches: 


SecleRibSpcO 

SecmcRibSpcO 

SecteRibSpc() 


leading edge ribs 
midchord ribs 
trailing edge ribs 


The ribs do not contribute to the 
calculation of structural weight 


wing's bending inertia, but are important for the 
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Plate-like Elements 


Core Type: 
Macro-core 
or 

Ti-core 


| ' i 


! I 


1 I 


Skin Thickness 

T f T 


Outer Thickness 

* 


Beam-like Elements 



S "ictaesf by three P™™** skin ihickness. 
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Bending Inertia 

The subroutine CompSec computes the structural moments of inertia and desired 
mass properties of each user-defined section. The routine moves sequentially through each 
section element, adding the contributions of each element to a set of running totals. Letting 
/ represent the section elements (e.g. leading edge skins, forward spar, etc.), and 5, be the 
area of the element, the routine keeps the following sums: 
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where £ is the Young's modulus of each area element, dS. From these totals, the program 
can extract the section's neutral axis location. 


SecCeruroidx = 


SecCeruroidz - 


centroid 


v centroid 


(ES), 

(ES) 

(£S), ’ 

(ES) 


and thus obtain the moments and products of inertia, 

(EI) x .-(ES\ l -z ctnln jES\ 

(EI) ZZ -(ES)„ -x cattn jES\ 

(£/) - ^ ES)^ z emln J ES) X ~ x cauro J < ES) Z + x cmm J- emt roii EE ) 
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Because (EI) a will generally be non-zero, the program must account for the inertia mixing 

between axes. We are interested in bending deflections about the x-axis, so we form the 
combination. 


SecEInertia = (El) - - la i 7 

The desired mass properties are calculated in a similar manner: 

SecLineDensiry = \ m = u p dS 

( X J, " 

(U-IJp 


p zdS 


SecCGx = x, 


CO 


SecCGz = z, 


(o. 

K 

(K), 


CG 


where SecLineDensiry is the mass of the structure per unit length of span and SecCGx and 
SecCGz represent the center of mass location. 

Bending Deflection 

Calculation of the wing bending moment and deflection is accomplished by 
analytical integration of the applied loads and wing properties. We start be defining the 

total load applied to the wing at any spanwise location as the wing weight per unit span 
minus the supporting force per unit span: 

P(y) “ Weight(y) - Support(y) 

We assume that this load varies linearly between defined stations, so that the load may be 
written as 
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Ky)-r!* + (y‘ 


y.) 


r.f - /*» 

Ay, ' 


fory t <y<y lH , 


with the index i running from 1 to NStations. We have assumed the following definitions: 


P* - Weight 1 * (y.) - Support 1 '*^) 
K ik> * Weight **^ ) - Support ^(y, ) 

Ay< -y, 


Now we can calculate the shear loading and bending moment applied to the wing by 
integrating the linearly interpolated P(y): 


Shear(y) - £ P{y')dy' 

Shear , * 0 

+ IK**** 

Shear { « Shear t + — *■ ^Ay, 


Moment (y) - P P(y’)y'dy 
Moment , - 0 

Moment ■ Moment, + Shear Ay + f ~ 

' ‘ V 2 6 


'Ay , 2 


It is physically reasonable to assume that the inverse of the wing bending inertia varies 
linearly between the stations where the wing section has explicitly been defined: 


(£/)(y) - 



Ay, 


T l 

■J- 


fory. <y<y. n . 


The second derivative of the wing deflection is just the ratio of the bending moment to the 
bending stiffness, or 


d 2 Pefl ( Moment {y ) 
ay 2 W “ El(y) 
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Combining the linearly intopolaled formula for the stiffness with the above moment 
definition, we obtain the following polynomial solution for the second derivative of the 
deflection: 


d 2 Defl \ 

< dy 2 ) x 


0 


f tf_Def£ \ 

< d y 2 L 


** + 2.(Ay,) 2 + + *jAy,) 4 


where the polynomial coefficients are given by 


*0, - Moment '({EI)?)' 1 

K t , - v a-fiar) 


„ l^f((£0,T) ~ ED *) 1 /**»'- p* , 

K u m -r- L L i *—!— |+ -^1 *J—t/ E r{+X' 

2 [ Ay ( J 6A y t (( )i > 

K m f-feggr 1 

6Ay f [ Ay, ] 

Now the first derivative and deflection itself can be found by analytic integration: 

’ °* ♦p§*) (*»)♦ *. ^ - *„ 

°y A 2 6 12 31 20 41 in 


The constants of integration (which determine values for the deflection and its first 
derivative at /=!) are chosen such that both wing tips are at z = 0. 
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Results 

Numerical results for the baseline airplane are shown in figures 5 through 8. 
Figure 5 shows the weight distribution and wing stiffness as a function of span. Note that 
the wing is heaviest in the region of the fuel tanks, outboard of the pressurized passenger 
cabin. Note, also that the weight distribution is slightly asymmetric. The only asymmetry 
in the baseline configuration (as it is described in the database) is location of the APU. 

Figure 6 shows the wing bending moment and deflection due to an elliptical lift 
distribution. Note in particular that the bending moment is negative at the wing root, due to 
the fact that the majority of the weight is carried outboard. 

Figure 7 shows the wing bending moment and deflection due to 1-g static loads 
applied at the landing gear. Because the airplane is nearly span-loaded during flight, the 
peak bending moment due to ground loads is a factor of 10 higher than the peak moment 
due to flight loads. Ground loads, however, apply their peak moment to the stiffest part of 
the wing-inboard m the pressurized passenger cabin, whereas flight loads apply their peak 
bending moments outboard where the wing is le ss stiff. 

Figure 8 illustrates the magnitude of these bending deflections, compared to the 
scale of the airplane. In each case the total wing-tip deflection is about equal to the 
thickness of the wing at the root 

Table 1 is a complete airplane weight breakdown. The structural weight is 
computed by the program; the weights and CG locations of all fixed equipment are fixed 

entries in the database. The program's estimate for the structural weight of the baseline 
airplane is 1 12,345 lbs. 
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—A— StalLinOens 
□ StarUnDens 



■A — StalElnartia 
O StarElnertia 


Figure 5 Wing weight distribution and bending stiffness. Note that the win® 

ofth^ianH thC locatl °" of fa eI tanks (between stations 4-5 and 1 1-12) outboLd 
of the landing gear, and outboard of the stiffest part of the wing. ' outtK>ar<1 
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StaMomant 


Span wise Index 


Span wise Index 


Bending moment and deflection for the wing subjected to 1-g flight loads The 





Bending 



feSh'rZ Ben <iing moment and deflection with the wing subjected to 1-e around 
loads. The wing is supported soley by its landing gear, located quTtefer inSS 
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Bending due to 1-g Flight Loads 
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